On the correlation energies for two interacting electrons in a 

parabolic quantum dot 



Omar Mustafa 

^-v . Department of Physics, Eastern Mediterranean University 

^^ i G. Magusa, North Cyprus, Mersin 10 - Turkey 

^^ ' email: omar.mustafa@emu.edu.tr 

^ ; PACS; 73.63.Kv, 73.63.-h, 73.90.+f 

(March 12, 2009) 
O 



O 

C/3 



C^ 



o 
o 



X 



Abstract 



the Coulomb coupling strength A relative to the confinement. 



C^ \ The correlation energies for two interacting electrons in a parabolic quantum 

dot are studied via a pseudo-perturbation recipe. It is shown that the central 
— ^ ■ spike term, (m^ — l/4)/r'^, plays a distinctive role in determining the spectral 

Ph ' properties of the above problem. The study is carried out for a wide range of 

o 

(N 
> 

o 



Advances in semiconductor technology have made it possible to fabricate ultrasmall 
structures that confine electrons on a scale comparable to their de Broglie wavelength. 
Structures that restrict the motion of electrons in all directions are called quantum dots 
(QDs). The simplest QD consists of an insulator, e.g. AlGaAs, and a semiconductor, e.g. 
GaAs. Between is a potential difference that confines injected electrons to a thin layer at 
the interface, in which electronic motion in the direction perpendicular to the sandwich is 
essentially frozen out. Effectively, QDs are considered two-dimensional systems demonstrat- 
ing typical quantum effects such as discrete energy levels and interference. As a realistic 
and computationally convenient approximation, a harmonic shape of the laterally confined 
potential ( i.e., a two-dimensional oscillator effective in the plane of the dot) is often used 

Several experimental [2-6] and theoretical [7-23] methods were invested for the study of 
spectroscopic structure of interacting electrons in a harmonic QD. For example, electron 
correlation has been investigated by the many-particle Schrodinger equation [20], by per- 
turbation [21], shifted 1/N expansion [17], and by WKB treatments and exact numerical 
solutions [16]. QD helium ground state in a magnetic field is obtained by a Hartree, Hartree- 
Fock, and exact treatments [22]. Magnetic field dependance of electron energies in QD is 
studied by the decoupled approximation [23], screening of an ionic potential in QD [24], ■ ■ • 
etc. 

In this paper we consider the simplest nontrivial problem of two interacting electrons, 
with effective mass m*, in a harmonic quantum dot. The Hamiltonian of which is known to 
decouple into an exactly solvable center-of-mass Hamiltonian 

Hr=^ + -Mul R\ (1) 

and a non-exactly ( or, at best, conditionally-exactly) solvable Hamiltonian 

Where, uJq is the characteristic frequency of the parabolic confinement, R = [fl + f2)/2, 
Pr = Pi + P2; M = 2171*, r = f\ — 7^2, Pr = {pi — P2)/2, and /i = m* /2. Moreover, when a 
magnetic field B is applied perpendicular to the plane of the dot one would simply amend 



the above Hamiltonians and replace Uo by the effective frequency Co = Juol + ool/A, where 
ujc = eB/m*c is the cyclotron frequency, and the spin energy term Eg = g*fiBBSz ( with 
Sz = [I — (— 1)™]/2, g* is the Lande factor, and hb is the Bohr magneton) could be added 
to the total energy of the dot. However, upon the substitution r = V2loq, Schrodinger 
equation for the relative motion Hamiltonian (2) eventually reads 



d^ m^-l/A 2 A 
dq'^ g^ q 



'^k,m{q) = Hfe,„^fc,m(g)- (3) 

Where S/j^ = 2Ek^m/{hu;o}, E^^rn is the eigenvalue of Hr in (2), A = \^lo/a*, a* = 



h?t/{m*e'^) is the effective Bohr radius, and lo = Jh/(m*uJo) is the characteristic length 
of the harmonic confinement. 

It is well known, on the other hand, that results from perturbation theory are limited 
to the case where A <^ 1 [10,21]. Moreover, the shifted 1/N expansion technique (SLNT) 
[17] and WKB treatments [16] lead to dubious accuracies in connection with level ordering 
and energy crossings [16,25]. However, results from exactly solvable potentials can be used 
in perturbation and pseudo-perturbation theories, or they can be combined with numerical 
calculations. Nevertheless, in the simplest case, analytical and semianalytical calculations 
can aid numerical studies in areas where numerical techniques might not be safely controlled 
[26-31]. 

In this work we shall use a pseudo-perturbation theory (PSLET) to study the spectro- 
scopic structure of the relative motion Hamiltonian of two interacting electrons in a parabolic 
quantum dot, eq.(3). Where, in the numerous methodical predecessors of a subset of papers 
[28-31], an alternative possibility has been sought in the power-law asymptotic expansions 
using some other small parameter. It has been noticed that the presence of the central 
spike, e.g. (m^ — l/4)/g^ in (3), just copies the effect of the centrifugal and/or centripetal 
force and immediately inspires the use of small shifted inverse angular momentum quantum 
number. An exhaustive description of the necessary formulae of PSLET accompanied by 
the persuasive verifications of their numerical usefulness ( by immediate comparisons of its 
results with available brute force numerical data) could be found in ref.s [28-31]. 

PSELT recipe starts with the augmentation of m"^ — 1/4, of the central spike in (3), by 
^z)(^D + !)• Where lo = I + {D — 3)/2, D is the dimensionality, / is angular momentum 
when D = 3 and I = \m\ when D = 2 (the dimensionality under consideration). Then, 
we simply use 1// as a pseudo-perturbation parameter, with I = l^ — f3 and /9 is a vital 
shift as it removes the poles that would emerge at / = for D = 3. Equation (3), with 
^(9) = (?^ + '^/Q')/2, therefore reads 



1 rf2 p ^ (2(3 + 1)1 + (3((3 + 1) l\,, A , ,, ^ . .X 



(4) 



where Ek,i = '^k,i/2 and Q is a constant that scales the potential V{q) at large - Id limit and 
is set equal to P at the end of the calculations, for any specific choice of Id and nodal zeros 
k. Next, we shift the origin of the coordinate system through x = l^^'^(q — qo)/Qo, where Qo 
is currently an arbitrary point to be determined below. Expansions about this point, a; = 
(i.e. q = go) would lead to 
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where 



v^^\x)=B2X^ + '^^^, 



(7) 



v'^''\x) = 5„+2 a;"+2 + (_i)- (2/^ + 1) 
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It is then convenient to expand Ek^i as 



oo 
n=-2 



(10) 



Equation (5), along with (6-9), is evidently the one - dimensional Schrodinger equation for 
a harmonic oscillator Vt^x'^ /2, with Vl? = 2B2, and the remaining terms in Eq.(5) are con- 
sidered as an infinite power series perturbation to the harmonic oscillator. One would then 
imply that 
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Where qo is chosen to minimize E]^i , i. e. 



dqo 
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and 
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Equation (13) in turn gives, with / = a/^. 



lD-P = ^qlV'{qo) 



(14) 



The shifting parameter (3 is determined by choosing lEf, ^ =0. Hence 
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where primes of V{qo) denote derivatives with respect to qo- Then equation (5) reduces to 
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Setting the wave functions with any number of nodes k as 



^fc,/(a;(g)) = Fk^x) exp{Uk,i{x)), 



(17) 



equation (16) readily transforms into the following Riccati equation: 
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(18) 



where the primes denote derivatives with respect to x. It is evident that this equation admits 
solution of the form 



uU^) = E urix) r-/^ + E Gkix) r("+^)/^ 



n=0 



n=0 
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Substituting equations (19) - (22) into equation (16) implies 
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The solution of equation (22) follows from the uniqueness of power series representation. 
Therefore, for a given k we equate the coefficients of the same powers of / and x, respectively. 
Although the energy series, equation (23), could appear divergent, or, at best, asymp- 
totic for small T, one can still calculate the eigenenergies to a very good accuracy by forming 
the sophisticated Pade approximation 
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(24) 



to the energy series (23). The energy series is calculated up to E^i /l^^ by 



Ekd = I E. 
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and with the P5^(l//) Pade approximant it becomes 



72 77(-2) 



EkAlO,9] = l'Eir> + Pi{l/l) 



(26) 



Following the above procedure, PSLET results are compared, in table 1, with the exact 
numerical ones ( obtained by direct numerical integrations, DNI) [16] for A = 1 and A = 10. 
To avoid exhaustive numbers of tables we do not list Garcia-Castelan et. al's results [16] 
from WKB, WKB single-parabola (WKB-SP), and WKB double-parabola (WKB-DP). In 
contrast with the WKB, WKB-SP, WKB-DP [16], and SLNT [13] resuhs the comparison 
between PSLET and DNI results implies excellent agreement. 
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In order to make remediable analysis on the effect of A, hence of the characteristic length 
lo (A ~ /o), we hst ( in tables II-IV) PSLET results for A; = 0, 1, 2 and A = 0, 1, 2, 4, 6, 8, 10, 12 
at different values of \m\. They are also plotted in figure 1. 

Figure 1 (along with tables II-IV) shows that the degeneracies associated with the har- 
monic oscillator confinement at A = are only partially lifted as A increases from zero ( of 
course, such degeneracies would completely be lifted when a magnetic field is applied per- 
pendicular to the plane of the dot). It also shows that the equidistance form of the energy 
levels at A = changes in the following manners; (i) for a given k, the spacing between two 
successive |?72| states decreases as A increases, and increases as \m\ increases for a given A, 
whilst (ii) for a given A, the spacing increases as the nodal quantum number k increases. 
One should nevertheless notice that (iii) s-states (with m = 0) shift up more rapidly than 
states with \m\ > 1, and for |?72| > 1 states with lower \m\ shift up faster than states with 
higher \m\ as A increases from zero. 

The above mentioned features (i)-(iii), in fact, build up the sought after scenario for the 
change in level ordering, that manifests energy crossings and spin-singlet [S^ = 0) spin- 
triplet {Sz = 1) oscillations, and inspires the vital role of the central spike term in (3). More 
specifically, the twofold nature of the central spike term in the effective potential, of eq.(3), 

T^ / X m^ - 1/4 , A ,^„, 

Veffiq)= ^,' +<f + - (27) 

explains the energy crossings as follows; (a) for m = it represents an attractive core that 
strengthens the confinement g^, whereas (b) for \m\ > 1 it represents a repulsive core which 
renders, along with the Coulomb repulsion, the potential less potent. This is why, for a 
given k, the energy of a lower \m\ state increases much faster ( more rapidly for m = 0) than 
that of a higher |m|, as A increases, and catches up with it ( hence energy crossings and 
singlet-triplet spin oscillations occur, or, at most, energy levels clustering is manifested). 
On the physical sides, the two electrons are farther apart for higher |m|. Moreover, for a 
given k energy crossings are not feasible between the corresponding states with different 
\m\. Whereas, states with a given k and |m| cross with states at lower k and higher |?72|. 
Therefore, the lowest three states (0,0), (0,1), and (0,2) never cross any other state ( i.e., 
they can never be depressed into a lower /c-state). 

The effect of correlation, between two interacting electrons in a harmonic QD, is therefore 
clear in the full energy spectrum for A > with all (A;,|m|)-states for the relative motion as 
shown in figure 1 and documented in tables I-IV. However, it should be noted that the level 
ordering reported by Garcia- Castelan et. al [15] is now changed, namely for the (0,4) and 
(1,1) states. Moreover, the (1,4) and (2,1) states seem to change order as A increases form 
12. 

To sum up, we have used a pseudo-perturbation recipe (PSLET) to study the character- 
istic length effect on the correlation energies for two interacting electrons in a parabolic QD. 
We have proved PSLET persuasive numerical reliability in comparison with direct numerical 
integration method ( in table 1) for A=l, and 10. Next, we have obtained the correlation 
energies for A = 2,4,6,8, 12 and k = 0, 1,2, and documented ( through figure I) that the 
level ordering reported by Garcia- Castelan et. al [15] is not absolute but continually bound 
to change as A increases from zero. Finally, the almost forgotten twofold effect of the central 

7 



spike term, in the effective two-dimensional potential (27), is now clarified to inherit a major 
responsibility for energy crossings. 
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TABLES 

TABLE L Comparison of PSLET energies (in fwOo/2 units) and the exact ones from direct 
numerical integration [16] for A=l and 10. 





A=l 






A=10 




(fe, \m 


) Exact 


PSLET {h 


:, \m 


) Exact 


PSLET 


(1,7^ 


20.3587 


20.3587 


fl 7' 


23.5040 


23.5040 


(0,9] 


20.3280 


20.3280 


(0,9] 


23.2188 


23.2188 


(3,2] 


18.5351 


18.5351 


(3,2] 


23.0339 


23.0339 


(2,4] 


18.4388 


18.4388 


(2,4] 


22.2217 


22.2217 


(1,6] 


18.3843 


18.3843 


(3,1] 


21.8721 


21.8715 


(0,8] 


18.3472 


18.3472 


(1,6] 


21.7355 


21.7355 


(3,1] 


16.6498 


16.6498 


(0,8] 


21.3954 


21.3954 


(2,3] 


16.4895 


16.4895 


(3,0] 


21.3513 


21.3140 


(1,5] 


16.4163 


16.4163 


(2,3] 


20.6504 


20.6504 


(0,7] 


16.3701 


16.3701 


(1,5] 


20.0186 


20.0186 


(3,0] 


14.9850 


14.9881 


(0,7] 


19.6037 


19.6037 


(2,2] 


14.5646 


14.5646 


(2,2] 


19.2438 


19.2438 


(1,4] 


14.4579 


14.4579 


(1,4] 


18.3753 


18.3753 


(0,6] 


14.3983 


14.3983 


(2,1] 


18.1420 


18.1420 


(2,1] 


12.6961 


12.6961 


(0,6] 


17.8543 


17.8543 


(1,3] 


12.5154 


12.5154 


(2,0] 


17.6671 


17.6660 


(0,5] 


12.4340 


12.4340 


(1,3] 


16.8431 


16.8431 


(2,0] 


11.0848 


11.0883 


(0,5] 


16.1628 


16.1628 


(1,2] 


10.6024 


10.6024 


(1,2] 


15.4916 


15.4916 


(0,4] 


10.4814 


10.4814 


(0,4] 


14.5546 


14.5547 


(1,1] 


8.7594 


8.7594 


(1,1] 


14.4622 


14.4622 


(0,3] 


8.5485 


8.5485 


(1,0] 


14.0379 


14.0381 


(1,0] 


7.2340 


7.2362 


(0,3] 


13.0720 


13.0720 


(0,2] 


6.6538 


6.6538 


(0,2] 


11.7903 


11.7903 


(0,1] 


4.8553 


4.8553 


(0,1] 


10.8495 


10.8496 


(0,0] 


3.4952 


3.4968 


(0,0] 


10.4816 


10.4816 
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TABLE II. PSLET correlation energies (in nuJo/2 units) for k={), |m| =0,1,2,3,4,5 and 
A=0,l,2,4,6,8,10,12 



\m\ 


A=0 


A=l 


A=2 


A=4 





2 


3.4968 


4.6391 


6.4428 


1 


4 


4.8553 


5.6557 


7.1251 


2 


6 


6.6538 


7.2872 


8.4994 


3 


8 


8.5485 


9.0864 


10.1331 


4 


10 


10.4814 


10.9564 


11.8885 


5 


12 


12.4340 


12.8638 


13.7112 




A=6 


A=8 


A=10 


A=12 





7.9373 


9.2644 


10.4816 


11.6184 


1 


8.4599 


9.6938 


10.8496 


11.9425 


2 


9.6480 


10.7425 


11.7903 


12.7975 


3 


11.1440 


12.1226 


13.0720 


13.9947 


4 


12.7978 


13.6861 


14.5547 


15.4049 


5 


14.5429 


15.3599 


16.1628 


16.9525 
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TABLE III. Same as table 2 for k=l. 



\m\ 


A=0 


A=l 


A=2 


A=4 





6 


7.2362 


8.2945 


10.0462 


1 


8 


8.7594 


9.4879 


10.8608 


2 


10 


10.6024 


11.1913 


12.3314 


3 


12 


12.5154 


13.0233 


14.0173 


4 


14 


14.4579 


14.9110 


15.8031 




A=6 


A=8 


A=10 


A=12 





11.5189 


12.8317 


14.0381 


15.1665 


1 


12.1368 


13.3327 


14.4622 


15.5362 


2 


13.4252 


14.4772 


15.4916 


16.4721 


3 


14.9840 


15.9254 


16.8431 


17.7388 


4 


16.6773 


17.5343 


18.3753 


19.2009 
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TABLE IV. Same as table 2 for k=2. 



\m\ 


A=0 


A=l 


A=2 


A=4 





10 


11.0883 


12.0757 


13.7327 


1 


12 


12.6961 


13.3720 


14.6650 


2 


14 


14.5646 


15.1195 


16.2015 




A=6 


A=8 


A=10 


A=12 





15.1801 


16.4736 


17.6660 


18.7833 


1 


15.8856 


17.0418 


18.1420 


19.1937 


2 


17.2480 


18.2613 


19.2438 


20.1978 
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Figures captions 

Fig.l: PSLET correlation energies for two interacting electrons in a harmonic quantum 
dot vs A = \/2lo/a* (the ratio of the oscillator length l^ and the effective Bohr radius a* 
. The energies are normalized with the oscillator energy huJo/2. The full lines represents 
states with /c = 0, dashed lines for /c = 1, and dashed dotted lines for k = 2. 
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